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Abstract 

We rederive the C(a 5 ) shift of the hydrogen levels in the non-recoil limit (m e /mp — > 
0) using Nonrelativistic QED (NRQED), an effective field theory developed by Caswell 
and Lepage (Phys. Lett.l67B, 437 (1986)). Our result contains the Lamb shift as a 
special case. Our calculation is far simpler than traditional approaches and has the 
advantage of being systematic. It also clearly illustrates the need to renormalize (or 
"match" ) the coefficients of the effective theory beyond tree level. 
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1 Introduction 



The Lamb shift, the shift between the hydrogen 2S\/2 and 1P\i% states, is without any 
doubt the most well-known bound state application of radiative corrections. It has the form 
m e a 5 (In a + finite) where the finite piece contains the Bethe logarithm, a state dependent 
term that must be evaluated numerically. The log term is relatively easy to extract and its 
calculation is presented in many quantum mechanics textbooks. The finite contribution 
is much more difficult to evaluate because it requires the application of QED to a bound 
state. In this paper, we rederive the complete 0(a 5 ) corrections in the non-recoil limitQ 
(m e /mp = 0) using NRQED, an effective field theory developed by Caswell and Lepage 
0, as extended by Labelle to study retardation effects. 

To construct NRQED, one writes down the most general Lagrangian consistent with 
the low energy symmetries of QED such as parity and gauge invariance, etc. The first few 
terms of this Lagrangian are given by (we follow the notation of [0]) 



(tr .„ D 2 D 4 . qcr ■ B . q(D • E — E • D) 

-2-Fermi = T\lD t + - h Or—.. 7 + Op — h 6 D — — 

2m Sm d 2m 8m 2 

. iqcr ■ (D x E — E x D) 

+5 S- " o 2 " + • • 

8m 2 

+ same terms with ip e — > ip p + photon terms, (1) 

where ip e and tp p are two component fields associated with the electron and proton, respec- 
tively. There are of course many other interactions, including four-fermion interactions, 
however, as we will discuss below, those are not relevant at 0(a 5 ). For the photon, we use 
the Coulomb gauge which is the most efficient gauge to study nonrelativistic bound states 
since it permits to isolate the Coulomb interaction (which must be treated nonperturba- 
tively) from all other interactions (which can be treated as perturbations). In that gauge, 
the photon terms in Eq.(|TJ) are 

- ~ A F, V F^ + 5 VP ^-A (k)— 2 A°(k) - 8 VP ^-A\k)— 2 A\k) (<% - + . . . (2) 
4 157T m 2 157T m 2 k 2 



Using D = z(p — qA) and D t = d t + iqA , the NRQED hamiltonian can be written as 



f- + gA -/j-^-(p' + P ).A+^A.A 
2m Sm 6 2m 2m 



-5 F ^-a ■ (k x A) - 6 D -^- 2 k"A° + 6 s -^<r ■ (p' x p)A° 

2m 8m z Am z 

+S s ^- 2 k cT ■ (p' + p) x A - 6 S -^<T ■ (k! x A(^))A (fc 2 ) + . 
Sm z Am z 



1>, 
(3) 



1 In the literature, these corrections are also referred to as the Lamb shift. We will also adopt this 
notation in the rest of the paper. 



1 



and the photon hamiltonian can be written as 



n photon = 7-(E 2 + B 2 ) - 5 v *A°(k)— 2 A°(k) + Sy^A^k^A^k)^ - + . . . 
2 157T m z 157T mr k z 

(4) 

The Feynman rules of Eqs.(|3|) and @ are given in Fig. [1]. 

The photon propagator requires a special treatment. We use time-ordered perturbation 
theory (in which all particles are on-shell but energy is not conserved at the vertices) and, 
as mentioned above, we work in the Coulomb gauge. In a time-ordered diagram, Coulomb 
photons propagate instantaneously (i. e. along vertical lines in our diagrams since we choose 
the time axis to point to the right) and transverse photon propagate in the time direction. 
The corresponding Feynman rules are given in Figs.0 and ||. 

As is shown in |2[], the power of NRQED is enhanced if one separates "soft" transverse 
photons (having energies of order 7 = Z\iol ) from the "ultra-oft" photons (with E ~ 7 2 /^) 
because the counting rules differ for the two types of photons. In that reference, it is also 
shown that the multipole expansion can be applied to the vertices containing ultra-soft 
photons. The distinction between soft and ultra-soft photons is particularly important in 
this work because the Lamb shift involves both types of contributions. We have there- 
fore separated the general photon propagator of Fig. [Q into a contribution corresponding 
to a soft photon, Fig. Pj(a)], and the contribution corresponding to an ultra-soft photon, 
Fig.[(|b)]. 

Finally, in NRQED loop diagrams, all internal momenta must be integrated over, with 
a measure d 3 p/(2n) 3 . 

The unknown coefficients 5r, 5p . . . can be found by imposing that NRQED is equiva- 
lent to QED for nonrelativistic scattering processes, i.e. by imposing 



NRQED scattering QED scattering amplitudes 

amplitudes expanded in powers of p/m. 

This is the so-called matching procedure. Notice that no bound state physics enters at this 
stage of the calculation. 

The coefficients appearing in (^) can be fixed by considering the scattering of an electron 
off an external field A^. This is illustrated, at tree level, in Fig.||. We have chosen the 
normalization of the interactions of Eq.(||]) in such a way that the tree level matching gives 
5=1 for all the coefficients appearing in (^). The first non-zero contribution to Syp comes 
from the one-loop QED vacuum polarization diagram, as illustrated in Fig. |] ] and, again, 
our normalization is such that 8yp = 1, to one loop J7|. We will still refer to this as "tree 
level matching" because only tree level NRQED diagrams are involved (similarly, n-loops 
matching will refer to the number of loops "n" the NRQED diagrams). The one-loop 
matching will bring 0(a) corrections to the coefficients 5's. 

We are now in position to proceed with the calculation. All NRQED calculations can 
be divided into three steps. The first one consists in using the counting rules to identify the 
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diagrams which will contribute to the order of interest. This first step not only permits to 
identify the relevant diagrams but it also fixes the order (in the number of loops) at which 
the coefficients must be matched. The second step consists in matching the coefficients to 
the order required and the last step corresponds to finally evaluating the NRQED bound 
state diagrams. 

In our case, we first need to isolate the NRQED diagrams contributing to order a 5 
in the non-recoil limit i.e. we neglect corrections suppressed by powers of m e /m p . As a 
simpler example, we will first isolate the diagrams contributing to order a 4 in the non-recoil 
limit (the full NRQED calculation of the 0(a 4 ) energy shift for arbitrary masses will be 
presented in Ref.||). In that limit, the only relevant diagrams are shown in Fig.Q. This 
can easily be checked by using the NRQED counting rules derived in Ref.[2]. Since soft 
and ultra-soft photons obey different counting rules, we consider their contribution in turn. 
Soft photons contribute to order : 

« ^Vc/, (5) 
m^rrip nip 

where p is the reduced mass and p and k are, respectively, the total number of inverse 
powers of m p and m e appearing in the NRQED vertices. In (|5p, we approximated p ~ m e 
which is exact in the non- recoil limit. The coefficient ( is defined by 

( = l + K + p-N TO P + Ylnh (6) 

i 

where Ntop is the number of intermediate state time-ordered propagators (see Ref. [2] for 
more details) and last term is the sum of factors of a contained in the vertices (including 
possible factors coming from the 8's). Now, to obtain the correction of order m e a 4 , we 
need to have p = and ( = 4 (see Eq.(|5])). The only way to have p = is to either have 
a Coulomb interaction on the nucleus line or no interaction at all. This already reduces 
greatly the possible diagrams. Turning now to the condition ( = 4, we obtain from Eq.(^]) 

k-aw + E^ = 3 - ( 7 ) 

i 

Since in first order of perturbation theory N top = 0, we are left with the condition 
K + Hi n i — 3- One possibility is k = 3 and n i — which corresponds to the relativistic 
kinetic energy vertex on the electron line. Another possibility is k = 2 and a i — 1 which 
can be fulfilled with the Coulomb vertex on the proton line and either the Darwin or the 
Spin-Orbit interaction on the electron line. There are no Coulomb interaction with only 
one inverse mass so the condition k = 1 and J2 n i — 2 cannot be satisfied. The three 
possible diagrams are illustrated in Fig.jT]. 

To order a 5 , and still in first order of perturbation theory, we must either increase K or 
TLi by one. It is not possible to increase the number of inverse electron masses k by one, 
but there are two ways to increase Yj n % by one - One is to include the one-loop corrections 
to the coefficients 5's of the vertices in Fig.0, so we will have to match these interactions 
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to one loop. The other possibility for k = 2 and J2 n i — 2 is to consider the new interaction 
corresponding to the vacuum polarization correction to the Coulomb propagator which is 
depicted in Fig.||. 

We now turn to diagrams in second order of perturbation theory, in which case Ntop — 
1. It can easily be verified that (, Eq.@, cannot then be made equal to 5. We have now 
uncovered all the diagrams containing only soft photons which contribute to order a 5 . The 
only remaining possibility is to consider diagrams with ultra-soft photons which lead to a 
contribution of the form given by Eq.(|5|) but now with (see Ref.[2]): 

C = Y / n l + l + p + K-N top + 2N^ + Y / M i , (8) 



where iV 7 is the number of ultra-soft photons in the diagram and Aii is the order, in the 
multipole expansion, of the i t h vertex and the sum is over the vertices connected to ultra- 
soft photons only. As before, we set ( = 5 and p = to obtain non-recoil corrections of 
order a 5 . For diagrams containing ultra-soft photons, the minimum value of Ntop is 1, 
because these photons propagate in the time direction (we again refer the reader to Ref. [2] 
for more details). Working at the zeroth order of the multipole expansion (Mi = 0) and 
considering only one ultra-soft photon (A 7 = 1), we then have the condition k + J2 n i — 3. 
Since the ultra-soft photon is necessarily transverse and transverse vertices contain at least 
one power of inverse mass (see Fig. [1]), k is bigger or equal to 2. In addition, J2 n i is at 
least equal to one (i.e. there at least a total number of one factor of a in the vertices). 
We therefore already fulfill the condition k + n i — 3 with the simplest diagram which 
corresponds to an ultra-soft photon connected to two p • A vertices, as represented in 

Pig-H- 

We have now identified all diagrams contributing to the order of interest. We now turn 
to the matching. From the above discussion, we see that we need to match to one loop the 
coefficients of the interactions contributing to order a 4 . 

To make NRQED agree with QED at the one loop order, we impose the relation illus- 
trated in Fig. ||10||. This matching was performed in JF| but, even though our final result is of 
course the same, our derivation differs sufficiently to be presented it here. The QED scat- 
tering amplitude (the left hand side of Fig. [0) can be expressed in terms of the usual form 
factors F\(Q 2 ) and F 2 (Q 2 ) in the following way (to be consistent, we use a nonrelativistic 
normalization for the Dirac spinors): 



u(p ',s') 
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F 2 (Q 2 )^ 



lo A°(Q)F 1 (Q 2 ) - — <r°M°(Q) Q J F 2 (Q 2 ) 
2m 



5m z 



Am 2 



JP ' x P)A° + 



2 AO 



Am 2 



Q M 



le 
2^2 



cr ■ 



(p ' x p)A° + . . . 



£ 5 



u(P, s) 



(9) 



where Q = p' — p, s and s' are the initial and final spin of the electron, respectively, and 
are the corresponding initial and final two component Pauli spinors normalized to 
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unity. We choose the convention that e is positive so that the charge of the electron is — e. 
In (|9]), as in the rest of the paper, we use m to represent the mass of the electron. Notice 
that the matching involves a double expansion. One expansion is in the coupling constant 
a and the other is the nonrelativistic expansion in Q/m which leads to renormalization of 
different NRQED operators. The nonrelativistic expansions of the form factors are |7fl 



F!(Q 2 ) 

F 2 (Q 2 ) 



3vr 



Q 



1 i m - 

ln( v 



+ . . . 



a Q 2 
a f , — — - + ... 



7T 



12m 2 



(10) 



where a e is the electron anomalous magnetic moment which, to the order of interest, can 
be taken to be a/(2ii). In the above Eqs., since Q° is of order of v 2 and |Q| is of order v, 
we have ignored Q° respect to Q. 

If we substitute (|l0]) in we obtain 



8a/, ,m. 3 



2a e + ... 



te 



4m 2 



(V.(p'xp)A»( 



l + 2a e + ... 



+ 0(Q 4 /m 4 ). (11) 



We must now compute the right-hand side of Fig. [I(J to complete the calculation of the 
one-loop renormalized NRQED coefficients. Since we are dealing with ultra-soft photons 
in Figs. 0(h)], fi~0|(i)l and fllQ"l(k)1, we use the special Feynman rules derived in Ref.[2]. 
Working at the zeroth order of the multipole expansion, Fig.[T(](h)] corresponds to 
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2epi 2ep'j, 
2m 2m ' 
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Oi ~ k 2 +A 2 
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2l 

2m 



2tt 2 
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k 2 + A 2 2Vk 2 + A 2 k 2 + A 2 



8a 
3^ 



In(^) - I 
A 6 



CeA i 



P ■ P 

4m 2 



(12) 



where, in the second line, we have used the fact that the integral is already proportional to 
p ■ p' to approximate p 2 ~ p' 2 ~ in the propagators. The corrections to these expressions 
will lead to higher order operators. 

Notice that we are only working with scattering diagrams when performing the match- 
ing, no bound state physics enters this stage of the calculation. To compute the amplitudes 
in Figs. [|i~0"l (i)1 and [|K](k)], we just need to evaluate the first one and then obtain the second 
one by replacing p by p '. For Fig.||T0"|(i)1, we can write 



d 3 k 



Sij - 



k 2 +A 2 



1 



2epj 2epj , 
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(13) 
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Here, E represents the on-shell energy p 2 /2m. Of course, in that limit, the propagator 
1/(E — p 2 /2m) is divergent, which signals the need for a mass renormalization. In NRQED, 
we perform mass renormalization exactly as in QED, i.e. we start by keeping E ^ p 2 /2m 
and subtract the mass counter-term: 



2 mi 



dkk 2 2k 2 + 3A 3 



2tt 2 



k 2 + A 2 2Vk 2 + A 2 
1 1 



Vk 2 + A 2 / E 



pi 

2//i 



(14) 



2m' 



— we get a series which, by construc- 



Expanding the term in the parenthesis around E 

tion, starts at order (E — p 2 /2m) 1 , which cancels the propagator l/(E — p 2 /2m) in Eq. (|13|) . 
One can then finally take the limit E — > p 2 /2m with for result, for the sum of Figs.|]T0|(i)l 
and 0(k)], 
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37rm 2 



(P + P 



rffck 2 
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5m' 



(15) 



Putting everything together, the complete right-hand side of Fig. |I(J is equal to the sum 
of Eqs.(0) and dH): 



5m 
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3tt 
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Now, by equating (]Tlf) and (0), we can evaluate 5 D and 5$: 

+ 2a e + 0{a 2 ) 



6s 



a 8 

1 + -7T 

7T 3 



ln(— ) H 

V 2A ; 24 



l + 2a e + C(a 2 ). 



(16) 



;is) 



Notice that, since the relativistic kinetic vertex does not enter the matching, we also 
findQ 



5 R = l + 0(a 2 ). 



(19) 



We are now ready for the third and last step of the calculation, the computation of the 
bound state diagrams per se. Since only 8jj and 8s receive an O(a) correction, but not 

2 One might expect that 5r be equal to one to all orders, since this interaction comes from Taylor 
expanding the relativistic expression for the energy, but this might not be so if the regulator (as is the case 
in this calculation) breaks Lorentz invariance. 
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5r, only Figs. 0(b)] and 0(c)] are needed for the 0(a 5 ) calculation. Let us now start with 
Fig. 0(b)]. 

r -iZe 2 f d 3 p'd 3 p T . /<7-p'xp\ T . ,. 
s • L 

= 5s ^< — > 

m(Zaf ( 5(j - (I + 1/2)) S(j - (/ - 1/2)) N 

" * 5 Anm + 1/2) I (7TI) l )^- 5 ^ 



(20) 



where ^(p) is the Schrodinger wavefunction (including the electron spin 3 corresponding to 
the quantum numbers n,j and I. 
In Eq.(|UD we used 

d3 P' c i P >. r ' (Tpxp' . apxr' ^ 



(2tt) 3 p' 2 4vrr' 3 

For the diagram of Fig. [0(c)], we find 

d 3 p'd 3 p , ( e(p ' — p^ 2 



dVd \- iP >m P ) 



(2tt)< 

= fc^(0)l' = fc=^. (22) 

Using the counting rules, we also found that the diagram depicted in Fig.[§] would 
contribute to 0(ma 5 ). This diagram corresponds to the well-known Uehling potential and 
is found to be 

[d?pWp , ( 1 , -(p'-p) 4 a 1 \ , , 

= 7 72^r* (p \- e W^W 5vp 15^ w^v? Ze r {v) 

= -5vpZ*£- i,|*(0)|^-^^^^o. (23) 
ibis m l 157T rr 

We finally turn our attention to the only remaining diagram, which is represented in 
Fig. [10]. The corresponding integral is (as shown in Ref. [2], in zeroth order of the multipole 
expansion we set p ' = p on the vertices): 

f d 3 kd 3 p 2e Pi 2ep 3 1 1 k t kj 

v ' n 2m 



In the non-recoil limit, the spin of the proton completely decouples from the problem. 
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(ftp f (p 2 
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(2,)3J — 



In a bound state, beyond tree level one must include an infinite number of Coulomb lines 
in the intermediate state. This can easily be seen from the counting rules, Eq.@. Indeed, 
adding a Coulomb line will not change ( because this increases both N top and J2i n i by one, 
which has no overall effect. Because of this, one must use the Coulomb Green's function 
for the intermediate state. Using the bra and ket notation, Eq . ([24]) must then be replaced 
by the well-known expression: 



^y2f dkk <nWo P \n'><n'\v\n>^ 

3tt ^ J E„- E„/ -k ( J 



This part of the calculation is well known and is carried out in many texbooks (see for 
example Ref.||). The result is 

Aa(Zaf\ if/ = ° 

AE = m^ [ -^L{ (26) 

where < E n > is the Bethe logarithm which can be evaluated numerically ||. Now, using 
Eqs.([l7|) and ([ID, we add the a 5 contributions from Eqs.(|0]) and(H) to Eqs.(H) and (g§) 
to obtain 



ln^— + ^ if I = 



4a (Za) 



J 2<E n , > 1 30 

AE = m—^^- { (27) 

37T \ , Z 2 ma 2 , 3 / <S(j-(l+l/2)) 6(j-(l-l/2)) v . f , / n 

111 2<E n j> ^ %{2l+l)\ (l+l) I > 11 1 / U ' 

which is the well-known Lamb shift. 



2 Conclusion 

We have calculated the complete O (ma 5 ) non-recoil corrections to the hydrogen energy 
levels, also referred to as the Lamb shift. The superiority of NRQED over the traditional 
approaches is twofold. Firstly, the calculation of the bound state diagrams is greatly sim- 
plified because the use of an effective field theory permits to disentangle the contributions 
from low and high momenta and only QED scattering diagrams need to be evaluated. Sec- 
ondly, the NRQED calculation is systematic in the sense that simple counting rules can be 
used to isolate the diagrams contributing to a given order in a, which is not possible in 
traditional approaches. 
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Figure 1: NRQED Vertices 
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